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CHAPTER I 


INTPOEUCTION 

1.1 Mptivatiu? and Relevance of Thesis 

Ihe lACATE (Lower Atnosphere Oonfxssiticn and Temperature E)^>eri- 
nent) ndssion vas a hi^ altitude balloon platfonn test vMdi employed 
an infrared rekdiometer to sense vertical profiles of the ooncentrations 
of selected atm^pheric trace constituents and temperatures. The con- 
stituents wesre measured inverting infrared radiance profiles of the 
earth's horizon. The radiometer line-of-si^t was scanned vertically 
across the horizon at approximately 0.25°/seoond. The relative vertical 
positions of the data points making up the profile had to be determined 
to approximately 20 arc seconds. 

The balloon systan for eicocirplishing the missic»i is shown in ' 

Fig. 1.1-1. It consisted of: (a) a 39 million cubic feet (zero pres- 

sure) balloon, (b) a loed bar oontaining the balloon control equipment, 
(c) a package oontedning additioned balloon control electronics with 
gondola recovery parachute, and (d) a gondola oontaining the research 
equipment. The balloon was designed to lift the paQ^load to a float 
altitude of ^^^proocimately 150,000 feet. 

Instrumentation to determine the attitude of the balloon platform 
consisted of a magnetcmeter and 3 orthogonally oriented precision rate 
gyros. The three rate gyros were employed to obtain an accurate time 
history of the angular velocity ocnponents of the r^earch platform for 
sul!»eguent data reduction and attitude determination. 

The main problem in the LACATE esqierimait is to determine the in- 
stantaneous orientation (i.e. , the attitude) of the instrumentation 
platform with respect to a local verticed. Moreover, this orientation 
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nust be determined to an accuracy of 1°. Cnoe this is kncwi, the orioi- 
tation of the line-of-si^t of the radioneter can be determined siiK» 
its relative motion with respect to the platform is prescribed. 

1.2 Present Status of the Problem 

Stablizing the balloon research platform or {predicting its orien- 
tation is a major prcblem ^diich nust be solved in all balloon borne 
experiments requiring line-of -sight instrumentation. Feedback oontrol 
systems have been used to stabilize the balloon platform with respect 
to an inertud reference frame Stability is obtained by suspending 
the platform at its center of gravit^^ and enploying some oontrol system. 
Oontrol system instnmentation includes sensoirs (rate gyros, digital 
star trackers, etc.) and reaction \i4ieels for torquing the platform. 
Systems of this type, however, are usually extremely ocnplex and costly. 

An alternate approach to this problem is to allow the badloon plat- 
form to swing freely fron its suspension point and then enploy seme 
method to determii^ its attitude (orientation) . Ihe orientation par- 
aaneters for the platform are determined by fitting the results obtained 
from a mathematicad model (vdiich simulates the balloon system) to those 
resvlts obtained from the platform's sensors (i.e., gyroscopes) . 

Severed numericad parameter estimation methods have been developed 
to determine the attitude (cprientation) parameters. The {moblem is 
normally solved by enploying an optimization ptrooess vhich minimizes 
the error between predicted and known output results. In the case of 
balloon research platforms the optimization problem involves the mini- 
mization of the svsn of the squares of the differences between the angu- 
lar velocity oenpenents obtained from the rate gyrosoepes and those 
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predicted fron the roatheraticzd model . However, vdth this approadi, 
the problan of determining tt» optimal decision variables, (i.e., ini- 
tied condition parameters) can require oonsid^able ocnputer running 
time. 

A promising, new e^proach for determining attitude of balloon re- 
seeuxh platformB involves observer state space reconstruction. For 
totally observable systems, state estimators can be constructed. The 

state estimator, vhidi is driven by adl plant ii^^uts ^nd outputs, can be 

. (3) 

used to determine the system state . 

Observer systems are state estimators constructed such that the 

error in the estimated state deca^ to zero over a finite time interval. 

By subtracting the plant model from the observer model, the error model 

for the reocnstructed state can be determined. This model consists of 

a system of homogeneous, first order dlfferentiad equations. The eigen- 

vzdues of the resulting eigenvalue problem can be chc^en such that the 

error decays to zero in a aiall interval of time. In this case, then, 

the estimates response apparoaches the actual states exponentially. 

1.3 Object of Thesis 

The two main objectives of this thesis are given as follows: 

1. Develop an observer model for predicting the orlGntation of 
balloon borne research platfoams. 

2. Bnaploy this obeerver model in conjunction with actual data ob- 
tained from NASA'S LACATf: mission in order to determine the 
platform orientation as a function of time. 

In order to achieve the above objectives it will be n ecessary to 
first develop a ^n^al three dimaisional mathematical model for simu- 
lating the motion of the balloon platform. This will be dis'cussed next. 
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CHAFTER II 

DEVELOfMENT OF BAIUOCN SYST^ MmiEMATICAL MOCCL 
2,1 Iciealizationa of Systan 

Die general bedloai system to be studied in this report is shc^n in 
Pig. 1.1-1. Die actual motion of this system is very ocirplex and in- 
volves various types of oscillations including bounce (verticed oscil- 
lation) , pen&LLations (ixiplane motion) and s{dn (rotation) . In general, 
it is necessary to first idealize this system before developing the 
mathematical model. For pur p o s es of this study, the fbllcMing ixlealiza- 
tions vd.ll be made: 

1. Die mass of the balloon, subsystems and interconnecting sub- 
systems, will be "Ivnfied" at the locations shOfei in Fig. 2.1-1. 

2. Die balloon will be treated as an "equivalent" rigid body. 

3. Die altitude of the balloon static equilibriun position (float 
atlitude) will be aMutned to be a constant during the entire 
period of observation; i.e. , changes in this altitude due to 
loesei or changes in the properties of helium will be neg lec ted. 

4. Die interconnecting cables will be considered to be infl ex i b le. 
Die above idealizations were applied to the general balloon systan 

shewn in Fig. 1.1-1. Die resulting idealized systan is shown in Fig. 

2 . 1 - 1 . 

Diere are two alternate approaches v^di can be followed for puz> 
pnoco ot modeling the balloon system; these are siimorized below. 

1. Die mathematical model for the entire balloon system (Fig. 

2.1-1) can be developed. Die major disadvantage of this ap- 
proach is that it requires knowledge of the aerodynamic forces 
acting on the balloon itself. Moreover, with this model, a 
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large number of generalized ocxsrdinates are needed to specify 
the oonfiguration of the system. 

2. An alternate approach is to develop the mti^matical model for 
predicting the motion of subsystems one and two (Fig. 2.1-2) . 
niis model does not include the aerodynamic forces acting on 
the balloon; hcMever, it does require information on tt^ motion 
of the radar reflector support point 0 (Fig. 2.1-2). The 
advanta^ of this model is that the number of (Agrees of freedom 
is decreased, and the aerodynamic force effects are automaticedly 
included if tha moticai of the support point is known. This is 
the model vAiidi will be used for carrying out the re^arch dis- 
cussed in this thesis. 

2 .2 Generalized Coordinates 

The generalized coordinates for a givai system are those coordin- 
ates which are employed to specify the oonfiguration of the system at 
any instant of time. In any mediariiccd system there will be as many 
generalized coordinates as there are degrees of freedom. In the case of 
the idealized Ixirped system shown in Fig. 2.1-2, six generalized coordin- 
ates are required to specify the balloon configuration. These are ocm- 
prised of six Euler angles which specify the orientation of the tvro 
subsystems. The three translational coordinates located at the radar 
reflector support point 0 are not considered to be geieralized coordin- 
ates, since these are known (prescribed) from data obtained frcm the 
radar tracking installation. 

In gaieral, the Euler angles give the orientation of the bo<^ oo- 

I I f 

ordinate axes (x^ ) relative to a fixed coordinate system (X^ ) . A 
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series of three rotations.edxut the body axis is sufficient to edlcw the 
bocfy axes to attain any (tesired ori^tation. 

Severed sets of Euler angles eue possible for fixing the orientaticxi 
of the subsystems. One set of Euler angles (Set III) was employed in 
this woric and is shown in Fig. 2.2-1. The sequeioe of the three rota- 
tions which define this set is described belcw. 

a. a positive rotaticai about the axis resulting in the X^’ 
body system, 

b. a positive rotation 6 about the axis resulting in the XV 
body system, and 

c. a positive rotation ip about the X^' axis resulting in the X^ ' ' 
body system. 

The transfomation equation for the above sequence of rotations is 
given as follows; i.e. , 

_i I • 

X = AX 

'vhere 


A = 


(c (a^)c (jrf) -s (0) s (Oi) ) (c i^) s (jzJ) +s (e)c i^) s m ) (-s (Oi)c (e) ) 


(-C(0)S(0)) 


(c (0)C (0) ) 


(S(0)) 


(S (0»)c («)) +S (0) S (pOc (Oi) ) (S m S (»ti) -S (0)C (0)c (0<) ) (c (0)c (4;) ) 


X denote the fixed system axes, and 
X' ' ' denotes the fij^ body system axes. 


2.3 Lagrange's Equation 

The mathematical model for sinulating the motion of the research 
platform will be develc^sed by aiploying Lagreoige's equation. The genend 
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fom of this equation is gi\^ as follows; i.e. , 

4 — -|^ = Q. (i-1, ...,n) (2.3-1) 

dt 3^^ 3q^ 1 

%<here 

L = T-V = Lagrangian/ 

T = kinetic eansrgy of the system, 

V = potential energy of the system, 

= generalized coordinates, 

= generalized velocities, 
n = the nuRber of gsiercilized coordinates, and 
= the nonoonservative gsiercdized forces. 

For purposes of this ^rk the friction at the si^:port points 0 and 
1 vail be neglected. In addition, the aeroc^niainic drag forces acting on 
subsystems 1 and 2 will also be neglected. Hence, the generedized forces 
are equal to zero and Bq. 2.3-1 reduces to 

^ — n o ■j— 

at ^ 

2.4 Kinematics 

In order to obtain the kinetic energy of the system, it is first 
necess^uY to develop the kinematic expressions for the velocity (angular 
and linear) of the subsystems. By enplcying the Euler angles, the 
angular velocities for subsystem (1) are given as follows; i.e., 

S(e^) + ij^), and (2.4-1) 

fli 

^^3 “ ((z^iC( 03 ^)c(ai^) + s(a«j^). 
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vrfiere 

« ocnponent of the emgular velocity of subsystem (1) 
along the i th botty axis (i = ?.,2,3) , and 

angles of rotation for subsi^tem (1) . 

Similarly, the angular velocities for subsystem (2} az?e: 

= (-^2 S( 4 < 2 )C(e 2 ) + 62 C(»i»2))/ 

V»2^^^ * ®^®2^ + and (2.4-2) 

= (1ZI2 0(82)0(012) + «2 ^^*^2^' 


where 

= oarpor^t of the angular velocity of subsystem (2) 
along the i th boc^ axis (i = 1,2,3) , and 
02,O^2»J^2 ~ angles of rotation for subsystem (2). 


By enploying snail angle e^pproximations (i.e., S(8) = 8, C(8) = 1) , 
and by r^lecting second order terms in 8 and 'I'f 0?*(2.4-l) and (2.4-2) 
can be written as follows: 


and 


w. 


w. 


w. 


(1) = 
(1) = 
( 1 ) _ 


®1' 

^1' 


w. 


w. 


w. 


( 2 ) _ 


( 2 ) > 


82, 

^2' 


( 2 ) _ 


= 


(2.4-3) 


(2.4-4) 
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The tr 2 uislaticxial motion of the support point 0 (Fig. 2.1'*2) is 
referred to an aoda vhich is fixed in space. Ihe translational vel~ 
ocity expression for this point is given as follows; i.e., 


where 


( 0 ) 

( 0 ) _ 


Xi, 




(2.4-f)) 


( 0 ) 


X3, 


V 


( 0 ) _ 


= absolute velocity ocnponents of support point 0 

(i = 1,2,3) cQong the X^”' bc3dy axis of subsystem (2) 

The velocity e^qsiressions for support point (1) (Fig. 2.1-2) are 
giv>ai as follows; i.e. , 


Vi^^^ = (Xi - ri4>i), 

^2^^^ = (X 2 + ^ (2.4-6) 

= (X3), 

vAiere 

= absolute velocity oonponents of point (1) alc»^ the 
body axis of subsystem (2) , and 
r^ = distance between point (0) and point (1) . 


The velocity expressions for point (2) are given as foUore; i.e. , 


where 


(2^ • 

^1 ' ■ ^ 1 ^ ■ ^ 2 ^^^' 

V2^^^ = (Xj + (2.4-7) 



= absolute velocity ccnponents of point (2) alcx»g the 
body axis of subsystan (2) , aixl 
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■ distance between point (i-1) and point (i) . 

Hie above equations were developed by enploying small angle approx- 
imations (i.e., S(8) « 8r C(6) > 1 ) , and neglecting seocfid order terms 
in € and 4^. Moreover, it also assumes that the nature of the ring and 
clevis sipport at point (1) is siudi that the difference between the spin 
angles and ^2 anall. 


2.5 System Lagrangian 

Hie general kinetic expression for subsystems (1) and (2) is given 
as follcKfVs: 


y&iere 


-im. 


2 i 

. ^i3. (i)‘ 

+ ) 


(i-1,2) 




(2.5-1) 



kinetic energy of subsysten (i) , 


« mass of suli^ystem (i) , 

lil^ 1 ^ 2 ^ ^i 3 * nonents of inertia of subsystems (il along 



the (2) body axis, and 

ooiponents of the angular velocity of subsystem (i) 
along the (2) botfy axis. 


The total kinetic energy T of the balloon system is obtained by 
siitming Bq. (2.5-1) and substituting from Bqs. (2.4-3), (2.4-4), (2.4-6) 
and (2.4-7) ; this giv^: 

^ ^ ^( 1 ) ^ ^( 2 ) 

= + (^ 3 )^) 

• • • 5 • 2 ^23 • 2 


(2.5-2) 
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For purposes of developing Bq. (2.5-2) the nonents of inertia ^ 

I 22 r^lected and and m 2 were treated as point nasses. 

ihe system potential energy is due to the presence of the oonse]> 
vative gravitationail forces and is given aa follcws: 

V » (2.5-3) 

vhere 

= -nij^grj^C(0j^)C(<iij^) , and 

( 2 ) 

v' ^ = -m 2 g(rj^C(ej^)C( 4 ^j^) + ^2^(62 

The system Lagrangian (L ) , viiich is defined in Bq. (2.3-1) , is 
obtained by subtracting the total potential energy (Bq. (2.5-3)) from 
the total kinetic energy (Bq. (2.5-2)) and is given as follows: 

L = T - V 

= + (^+rj^e^)2 + (X^)^) 

+ (X2+rj^e^+r202)^ + (^3)^) + 

+ grj^C(0j^)c(4#j^) 

+ m 2 g(r3C(0j^)C(a<j^) + r2C(02)C(a^2))- (2.5-4) 

2.6 System Math Model 

The equations for the motion of the bedloon platform are obtained 
by sitetituting the lagrangian frcm Bq. (2.5-4) into Bq. (2.3-2). The 
resulting equations are given teXow. 
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m^lSl + 

( 2 . 6 - 1 ) 

" 21«1 * " 22«2 * •‘ 22®2 " ‘* 2 “ l ' 

’‘11*1 ■ ^1*2 ' 

(2.6-2) 

•• •« 

"'22*^2 ^22^*^ “ ^2^2 ' 

I 

•» 

- Of and 

(2.6-3) 

S2 “ 0^ 

where 

/ ^ 2 
"11 * ■*■ *"2^^1 ' 

"^2 * "^^2' 

"21 • ™2^2' 

T022 * 

^11 “ 

’‘22 ■ "’2 '^ 2 ' 

» (mj^-Hn2)r^f 

*2 - V2' 

12 ) 

= acceleration oortponent of point (0) along the Cj body 
axis, and 

( 2 ) 

^2 * acceleration ccnponent of point (0) along the e^ ' body 
axis. 

Bgs. (2.6-1) and (2.6-2) were developed by assundng anall displaoanents; 
i.e., c(e^) * c(o^^) « If s(©.) = e., s(0/^) « 4^^. 

Bq. (2.6-3) yields that is a constant. In this stucfy the pre- 
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cislon of this model was inproved by enplcying the transfomation equa- 
tion for the angular veloci^ . As stated earlier, it will be as- 
sumed that (2.6-3) will be r^lac^fi with the follow* 




(2.6-4) 


^ is the angular velocity obtained from the rate gyro 
( 2 ) 

mounted along the body axis. 

Bqs. (2.6-1) and (2.6-2) can be written in matrix fcxm as follows: 


MHi + Krii - Vi 


(2.6-5) 


Ei-oii 


(i=l,2) 


(2.6-6) 


4-2J ' 

^1 ”^2' 

"^1 ”22 


’‘u “ 


f 
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( 2 ) 


( 2 ) 

^ cxitfjotient of the angular velocity of subsystCDi 
( 2 ) 


(2) along the boc^ axis 


62 »w. 


( 2 ) 


W. 


( 2 ) 


2 - ^ ^ ooR|)Gnent of the angular veloci^ of subsystem 

( 2 ) 


(2) edong the 62 body axis 


C • (0 1), 




, arid 


e. 




Ihe nunerical valura for the and coefficients were ccxnputed 
by enploying the Enrpb-r'iies of the balloon system t^ch are given in 
Tzdale 4.1-1. these resulting values are given in Table 4.1-2. 

2.7 State Variable Form of Math Model 

Ek^. (2.6-5) and (2.6-6) can be expressed in general state variable 
form as follows: 

Lq«Nq + Ru, and (2.7-1) 


y • Cq 


(2.7-2) 
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where 



10 0 0 

0 10 0 

f 

0 0 lOj^ 

0 0 m^i m22 

~0 0 10 

0 oil 
-k^ 0 0 0 ' 

0 -^22 0 0 


R = control matrix = 



f 


u = control variable. 



( 2 ) 

W = the ocnponait of platform angular velocity along the 
system (2) bocfy axis, and 

C = (0 0 0 1) . 


Eqs. (2.7-1) and (2.7-2) can be fxorther expressed as follows: 


q = Pq + Bu, and 


(2.7-3) 



21 


v^nere 


y = C q 



(2.7-4) 


U»e elements of the A and B matrioes are givai belcw; i.e., 

_ ^2 ’"ll 
®31 3 ' 

k22 

^32 " 3 ' 

n»21 ^11 

Ni " 5 ^ 

_ "^1 ^22 
^42 " 5 ' 

^ "^2^1 " "*12^2 

b3 = j , 

^ _ "^1^1 "b.1^2 
^4 = 3 ^ 

^here 

3 = V "^2 ■ "bA2- 


The numerical vcdues for the and coefficients viere ocnfuted 
by employing the values of )c^j and viiich are given in Table 4.1-2. 
1hei» resulting VcLLu^ are givm in Tables 4.1-3 and 4.1-4. 



C31APTER III 


DEVELOPMEin' CF CBSEFKTEK SY^TTEM (OIHEMATICAL MXEI^S 


3»1 Oono^it: of Observability 

The cdsservability of a system inplies the deteminablLity of the 

system state fion an observation of the output over a finite time in- 

(4) 

tervcd starting from the instcmt at v^iich the state is desired . It 
is assumed that the system inputs, outputs and mathematical nodel are 
known. 

For purposes of refer^Kdng the work in this chc^Tter, the state 
variable form of the bcdloon's model (EJgs. (2.7-3) and 2.7-4)) will 
be rewritt^ below; i.e. , 


q = A q + B u, 
q (V = and 
y = Cq. 

v4iere 

q is the n th order state vector, 
is the unknown initial state vector, 
u is the single (scalar) input, 
y is the single (scalar) output, 
n is the order of the system, 

A is a n X n matrix, 

B is a n X 1 matrix, and 


(3.1-1) 


(3.1-2) 


C is a 1 X n matrix. 
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Gq. (3.1-1) determines the plant c^nendcs and indioatee hOM the 
input (or ccntrol) u affects the state vector q. The matrix A char- 
ac±erizes the plant <^nandcs vAien u is not pr^ent; this is the so- 
called free-respcnse case. The matrix B determuies how the plant re- 
sponse is effected by the ii^Jut (or oontrol) vector u. Bq. (3.1-2) 
defines the relaticnship between state vector q and the system output y. 

For system observability, the basic question is as follows: "Is 

it possible to identify the initial state by observing tiie output 
(y) over a finite time interval?" Precise definitions of system ob- 
servability are given as follows: 

1. Definition 3.1a. A state q^, i.e. q(t^) of a system is said 
to be observable at time t^, if kncxirledge of the input u(t) 
cBid output y(t) over a finite time t^ < t < t^, oonpletely de- 
termines the state c^. Otherwise, the state is said to be 

unobservable at t . 

o 

2. Definition 3.1b. if all system states q(t) are observable, 
then the system is said to be completely observable or just 
cbservable. 

3. Definition 3.1c. If the state q^^ is observable and if the 
knowledge of the input arx3 the output over an arbitrarily small 
interval of time suffices to determine q^ (ind^)endeit of t^) , 
then the state is said to be totally observable. 

4. Definition 3. Id. If all the states q(t) are totally observable, 

(5) 

then the system is said to be totally observable' / . 

The necessary oondition for observability of the balloon system is 
given in the following section. 
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3.2 Observability of Ballcxn Systenis 

The following theorem (Ref. 3) can be employed to determine the 
dsservability for general time - invariant systems. 

Theorem 3.2a. The time-invariant system described by Bqs. (3.1-1) 
and (3.1-2) is totadly observable if and only if the ooiposite matrix M 
has rank n, where 

C 

CA 

C(A^) 

M = 

— , and 

cuTh 

n, C and A are defined in Bqs. (3.1-1) and 3.1-2). 

The observability of the belloon system described by Bqs. (2.7-3) 

and (2.7-4) can now be determined by shewing that the cerposite matrix 

2 3 

M has rank 4. For this purpose, the matrices C, CA, c(A ) and C(A ) are 
given as follows: 

C = (0 0 0 1) , 

CA = (a^^ a^2 00)/ 

C(A^) =0 0 a^j^ ^42^' 

C(A ) = ( (a^j^a^^ + ^42^41^ ^^32^41 ^ ^42^ ^ * 

The specific form of the oonpMite matrix M for the balloon system 
model is given as follows, i.e. 





M 


CA 

C(A^) 

C(A^) 


; i.e. 


M = 


41 

0 




a 


42 

0 




a 


0 

0 

41 

0 


1 

0 

^42 

0 


In order to show that the rank of M = four, it is necessary to 
prove that the deterndnant is non zero. It can be shown that the dc ' 
termLnant for M is noi zero if the following expr^sion is non zero, 
i.e., if 


2 2 
Wl ^2 


0 . 


Thus, the balloon system is ooipletely observable since the composite 
matrix M has rank = 4. 


3.3 Observability of Balloon System With Output Bias 

Previous studies have been conducted for the lACATE system in order 
to determine the nature of the balloon's platform motion^®! The time 
history of the platform pendulation angles (6(t) and ii<(t)) was deter- 
mined by integrating the output from the rate gyros. The study indi- 
cated that the platform motion consisted of small oscillations siperim- 
posed on a line with (nearly) constant slope. These results sug^sted 
that the gyroscopes oont^dn a cxxistant bias error. 
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In order to take into account the error in the output (y) die to 
bias in the gyrosoc^ies, the matrices in Bgs. (3.1-1) and (3.1-2) are 
defined as follows; i.e.. 



( 2 ) 

w = the ocnporont of the angular velocity of the platform 
along the specified body axis, 

= oorrespoi:xiing gyroscope bias, and 

C = (0 0 0 1 1) . 

Uie observability of the system model described by EJqs. (3.1-1) 
and (3.1-2) can be ascertained by showii^ that the corresponding oom- 
posite matrix M has rank 5. It can be shewn that the corpe^ite matrix 
M for this system is given as foiled; i.e.. 



c 


M = 


CA 

C(A^) 

C(A^) 

C(A^) 


M = 


41 

0 


*42 

0 


^®41®3l‘^^42^41^ ^^41^32'^^42^ 


0 

0 

^41 

0 


1 

0 

*42 

0 


^®41^3l'^®42®41^ ^®41®32'^®42^ 


It can be shcMn that the detexminant of the Jbove matrix is non-zero 
if the resulting expression is non zero; i.e., if 

(r2+rj^) ^ 0. 

Thus, the system described by Bqs. (3.1-1) and (3.1-2) is oonpletely 
observable since the conposite matrix M has ran)c = 5. 

3.4 Develqpnent of Full Order Identity Observer for Balloon System 
Without Bias 

An n th order identity observer (or asynpotic-state estimator) can 
be constructed for the oorppletely observable n th order plant described 
by Bqs. (3.1-1) and (3.1-2). The observer is described by the follcwing 
equations; i.e. , 

z-FZ + Bu + Gy, and 

Z (t J = Z. 
o o 


(3.4-1) 



v^iere 


Z is an n order estimate of the state vector q, 
is an estimate of the unknown initial state vector 
G is an n X 1 matrix, 

F = (A - QC) , and 

Y, u, B and C are as described in Eqs. (3.1-1) and (3.1-2). 

An inspection of Bq. (3.4-1) reveals that the state estimators 
response (Z) vdll be detezmined from a consideration of the observer 
dynamics, external ii^xits and plant outputs, ihe observer dynamics are 
controlled by the F matrix. The external ii^t u oontribiites to the 
state estimator's response via the control matrix B. From Eq. (3.1-1), 
it can be seen that this control matrix B and the external inputs u are 
identical for both tte plant and observer. 

For cKxurate state space reoonstructian the plant outqwt y must be 
fed into the dDserver model. By coupling the plant out^xtt y to the ob- 
server via the G matrix, the observer loeoomes a closed loop estimator. 
This is illustrated in Fig. (3.4-1) . 

The F matrix in Eq. (3.4-1) is constructed sixh that the difference 
between the output of tl% observer model and the plant model is zero 
over seme finite time interval. This error (E) is given as follows; i.e. 

E Z - q. (3.4-2) 

Subtracting Bq. (3.1-1) from (3.4-1) yields the follcwir^: 

E = FE. (3.4-3) 

The solution of E)q. (3.4-3) yields: 
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Fig. (3.4-1) Block Dlagraun For Asymptotic State Estimator 
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E = E (3.4-4) 

o 

v4iere 

- %■ 

It is clear from Eq. (3.4-4) that the error, E, dece^ ejqxx^tiedly 
to zero if G is chosen such that all of the eigenvalues of F are nega- 
tive or have negative real parts. Also^ these eigaivalvies nust be more 
negative than the eigenvalues of A to insure accurate response. 

Ihe detailed form of the F matrix in Bq. (3.4-3) can l3e obtadned by 
substituting the form of the A and C matrices vhich are defined in Eq. 
(2.7-3) and (2.7-4). This yields the following; i.e. , 

F = (A — GQ , i.e . , 


vhere 
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Substitution of this into Bq. (3.4-3) yields the following eigenvalue 
problem: 


(F - XI)X » 0 


where 

I « the identxt^ matrix, emd 
F = matrix defined in Bq. (3.4-1) . 

Ihe necessary and sufficient ocndition for determining the eigen- 
vedxaes of the matrix F is that 


|; i.e, , 


-X 

0 

1 

^1 

0 

-X 

0 

(i-g 2 ) 

^31 


-X 

-^3 

^1 

^2 

0 

(- 94 -^) 


Ihe diaracteristic polyncmial cbtained by ejqianding Bq. (3.4-5) is 
given below; i.e. , 

(X^) + (g^X^) + ( (gia4i“a3i-a42'^2®42^ 

+ ( ^■g4^3i'^3®41^ ^"^2®41'‘'*42*31 

^2^2^4l"^2^2*31^ “ °* 

Critical damping of the error E is obtained by determining the G 
matrix sudi that all of the eigenvalues are negative and equal. This 
yields the following vzdues for the G matrix: i.e. , 
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93 * ((- 4 X^ 4 g^a 2 j^)/a^^) , and 

94 • (-4X) . 


3.5 Develoianent of Identity Otaserver Systera for Balloon System With 
Bias 

In the case of the fifth ozrder identity observer systan the matrices 
in Eq. (3.4-1) have the following form; i.e. , 



0 1 -9j^ -g^ 

0 0 (l-gj) ^2 

*32 ® -93 

*42 ° ■94 "94 

0 0 -gg -g^ 


^2 


G = 


^3 


i 


9 


4 
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and y, u, B and care as described in Eqs. (3.1-1) and (3.1-2). 

The eigenvalues for the F matrix corresponding to the bias model 
are obtadned fron the followij^ condition; i.e.# 
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-X 

0 

1 

-^1 

-^1 

0 

-X 

0 

(1-^2^ 

-^2 

*31 

*32 

-X 

■^3 

^3 

*41 

*41 

0 

(-X-g^) 

-^4 

0 

0 

0 

■^5 

(-X-gg) 


The characteristic polynomial of Eq. (3.5-1) is given as follows; 

+ ( (94*^5) + ( <9 j ^*4 i ‘^2*42‘*31“*42^ 

+ ( < 93 * 41 ”^ 4 * 31 “^ 5 * 31 "^ 5 * 42 ^ 

+ ( ^"^2*41*3l'^2*41*32‘^*42*3r®41*32^ 

^^V42"^l-^5*41^2^ “ °* <3.5-2) 

Ihe final fom 01 the G matrix elements (obtained from the oondi- 
tion for critical danping) are given as follows; 

g^ - <<^0^^“^2*42'*'*3l'^*42^/*41^ ' 

^2 “ <<5X^““42®3r®41^2>/<-*42*31-^*41^32)^' 
g^ * ( (— lOX <^3X^*42^ ^ ^*41^ > 

94 • (“5X -gg) , and 

g^ » ((-X )/«42*3l"*3i“*41^2^^* 
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CHftPIER IV 

RESOUTS AND OOWCIUSIONS 

4.1 13ata for lACATE Mission 

Figure (1.1-1) illustxates the actual lASSTE balloon system and 
Figure (2.1-2) illustrates the c o rresponding idealized system used in 
this stuc^. Ihe values for the various lengths and inass^ of the ide- 
alized system are given in T^le (4.1-1) . 

The ntnerical values for the elements of the M and K matrices of 
Bq. (2.6-5) were conputed based on the data given above. The resulting 
values are presented in Table (4.1-2) . The nunerical values for the A 
and B rnatrioes of Bq. (3.1-1) are giviai in Tables (4.1-3) and (4.1-4). 

The eigaivalue prcislem for the balloon system was ''olved analyti- 

2 

cally. The solution for the eigenvalues (Q ) j and corresponding eigen- 
vectors are presented in Table (4.1-5) . The vali:^ of Qj represait the 
natural frequencies of the system. The modal shape functions and per- 
iods corresponding to each natural freguenoy are shown in Table (4.1-6) . 

Results for the balloon observer system were obtained by enploying 
two s^aarate time intervals. The equations for predictino the boc^ axis 
aoceleratim ccnponents from sensor cxitputs are presented in ^^jpetidix A. 
Plots of the faotfy axis acceleration ocnponents over the two time inter- 
vals are given in Figures (4.1-1) through (4.1-6). 

4.2 Results From Simulation Study 

Bq. (3.1-1) was erployed to simul?te the balloon platfcm angxiLar 
velocity and angular displacement. The inputs and outputs employed with 
tte simulated system were of tlie same order of magnitude as predicted 
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TABLE 4.1-1 


Idealized LACATE System Prt^Jerties 

(distance from point 0 to mass nij^) = 75 ft. 

T 2 (distance from mass to m 2 ) = 15 ft. 

(linf)ed mass) = 135 Ib^ 
ni 2 (lunfjed mass) = 375 Ib^^ 


TABLE 4.1-2 


Coefficients of m and k Matrices 


= 89156 
m^2 = 13111 
= 13111 
m22 = 2622 
= 38278 
k22 = 5629.4 


(Ibf • • ft) 

(Ibf • s^ • ft) 

(Ibf • s^ • ft) 

(Ibf • • ft) 

(Ibf • ft) 

(Ibf • ft) 



-1.621 (s”^) 


^1 

= 1.1923 (s 


-2, 


‘32 


a.. = 


‘41 


a.„ = 


42 


8.107 (s"^) 
8.107 (s’^) 


TABLE 4.1-3 


Coefficients of A Matrix 


TABLE 4.1-4 


Coefficients of B Matrix 


.( 1 ) - 


-0.01334 (ft~^) 


( 1 ) - 


= 0.0 


= 

^3 

.( 2 ) 


0.01334 (ft’^) 


= 0.0 


36 





















ORIGINAL PAGZ IS 
OF POOR QUALITY 


















by the ac±ual lACATE flic^t data. Fbr conparison purposes, all of the 
sinuLation runs aiployed Identical initial oonditions aiui eigenvalues; 
the magnitude of the latter v;as equad to -0.6. 

The angular velocity (6) and angular displacement (6) predicted by 
the fourth order observer system are shown in Figs. (4.2-1) - (4.2-4) 
for the case when bias is not present in the output. The free response 
case is illustrated in Figs. (4.2-1) acnd (4.2-2), vdiile in Figures 
(4.2-3) and (4.2-4) , results are giviai for the case when a ramp irput 
was aiployed. 

Figures (4.2-5) - (4.2-13) illustrate the free response of the 
fourth order observer system for the case when bies is present in the 
output. The plant cxitput for the case shown in Figs. (4.2-5) and (4.2-6) 
contains a constant bias, vdiile a linear bias was used to obtain the re- 
sults shown in Figs. (4.2-7) and (4.2-8). 

Figures (4.2-9) - (4.2-13) give the free response results for the 
fourth order observer system vhen a high frequency bias >3 = 0.005 
cos(13t)) is present in the plant output. The results for the angular 
displacemeit from the fourth order observer system are presented in 
Figs. (4.2-9) and (4.2-10). Figs. (4.2-11) and (4.2-12) present the 
results for the angular velocity predicted by the fourth order observer 
model. Fig. (4.2-13) presents the results for the plant outpit (i.c. . 
y « + 6) . 

Figures (4.2-14) - (4.2-21) displey the free response results from 
tlv; fifth order observer system for the case v^ien bias in the output is 
pr^ent. In the case of Figs. (4.2-14) and (4.2-15), the output con- 
tained a constant bias, while a linear bias was used to obtain the re- 
sults shown in Figs. (4.2-16) and (4.2-17). 
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Figures (4.2-18) - (4.2-21) show the results from the fifth order 
observer system vAien the plant output contains a high frequency bias 
(qg = 0.005 cos (13t)). Pigs. (4.2-18) and (4.2-19) present resxilts 
for the angular displacement predicted by the fifth order dbserver sys- 
tem. Ihe angular velocity predicted by this system is shown in Figs. 
(4.2-20) and (4.2-21). The results for the plant output (i.e., y » 

+ q^) are presented in Fig. (4.2-13). 
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Fig. (4.2-5) Simulation Results For Angular Displacement 
(4th Order Observer, u=0, A=-0.6, Pias=1.0 x 
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4.3 Discussion of Siimlation Stu(^ Results 

The r^ults in Figs. (4. 2-1) -{4. 2-4) verify that, after a finite 
time interval , the fourth order observer system will reconstruct the 
plant state exactly if all plant inputs and outputs are known, ihe 
use of repeated eigenvadues (i.e. , ~ ~ 0.6) resulted in an error free 

re»X3nse after a period of 20 seconds. 

The fourth order observer results in Figs. (4. 2-5) -(4. 2-13) were 
obtained for the case v^en the plant output contained a bias (3(t)) . 

From Figs. (4.2-5)- (4.2-13) , it is seen that the reconstructed states 
contained an error resulting from the bias present in the plant output. 

TVro methods can be used to determine the steat^ state error in this 
reconstructed state due to cutput bias: 

1. Revise Bq. (3.4-3) to include the error due to bias and then 
solve this equation. 

2. Cdtpare the results for the reconstructed state to the au:±ual 
state directly by enploying the curves in Figs. (4.2-5)-(4.2-13) . 

For purposes of this work the latter method will be employed. 

For coiparing the reccnstructed state to the actual state, a relative 


error, ER^(t^) will be used. Ihis is defined as follow; i.e.. 


ER,(t.) = 


E?V^(tj) 


qilt^ 


(4.3-1) 


vhere 


EAi(tj) = magnitiKie of the error in the state at t = t^, and 
q^(tj) = maxirmm value of the actual state variable at t = t^. 
Ihe curves in Figs. (4.2-5) and (4.2-6) indicate that the error in 
the reconstructed state is a constant when the plant output contains 
a constant bias (i.e., S = 0.001) . In particxilar, for this case, a 5% 



relative error (ER ) in the outjHit resulted in errors (ER 2 and ER^) of 
17% and 100% in 0^ and 62 respectively. 

When a linear bias (i.e., 6 = 4.0xl0"^t) is present in the plant 
output, a linear error resulted in the reconstructed state (see Figs. 
(4.2-7) ard (4.2-8)) . At t = 43s, ER^ = 6 %, ER 2 = 23% ani ER^ = 133%. 

Ihe res\alts shown in Figs. (4.2-9)- (4.2-13) indicate that viien the 
plant output contained a high frequency bias (i.e., 6 = 0.005 cos 
(13t)), the angular displacement ( 6 ) and the angular velocity (0) pre- 
dicted by the fourth order observer contained a sinusoidcd error. Ihe 
error in the reccaistructed state at peak anplitudes are ER^ = 167%, 

ER 2 = 9% and ER^ = 33%. 

The results obtained from the fifth order observer are shown in 
Figs. (4. 2-14)- (4.4-21) . Figs. (4. 2-14) -(4. 2-16) indicate that, for the 
case when the plant output contains a constant bias, the actual state will 
be reconstructed exactly after a finite interval of time. 

Figs. (4.2-16) and (4.2-17) show that vAien the plant output contained 
a linear bias (i.e. , 6 = - 1 x 10 ^t) , a constant error resulted in the 
reconstructed state. At t = 43s., the resulting errors are given as 
follcws; i.e., ER^ = 153%, ER 2 = 64% cind ER^ = 100%. 

Figs. (4. 2-18) -(4. 2-21) show that when the plant output contcdns a 
high frequency bias (i.e., B = 0.005 cos (13t)) the state predicted by 
the fifth order observer contains a sinusoidal error. Hie errors in the 
reconstructed state at peak anplitudes are ERy = 167%, ER 2 = 6 % and 

ER. = 35%. 

4 

In genered, the results of the simulation study show that the form 
of the error in the reconstructed state will d^iend both on the order 
of the observer model and the type of bias (i.e., constant, linear, etc.) 
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present in the plemt output. When an nth order bias is present in the 
plant output, the fifth order observer will yield a more accurate response 
than the fourth order <±>s€rver. For this case, the error in the fifth 
order cA>server will be of order n-1, v^iile the error in the fourth order 
observer will be of order n. 

For a specific form of bieis, the magnitude of the errors in the re- 
constructed state will depend on the elements of the F and G matrices. 

For the particular ceise vdien = - 0.6, the magnitude of error in the 
reconstructed state Z 2 was much smaller than that in Z^. Moreover, for 
the case v^en the bias was of a sinusoidad nature (6 = 0.005 cos (13t) ) , 
the error in the reconstructed state Z 2 was less for the fifth order 
observer model; however, a more accurate response was obtained in state 
Z. for the fcxarth order observer. 


4,4 Ot«erver Results for Det 


Orientation of Balloon Platform 


Ihe fourth and fifth order observer system models were employed 
to determine the orientation (© 2 ) of the balloon platform in the 
plane for the time interval t = Os dnitial lACAlE data recording time) 
to t * 500s. The c^erver transient error was assumed to decay to zero 
after an elapsed time period of 250s. At this time the initial condition 
for integrating the gyroscope was set equal to the angular displacement 


(Z 2 ) predicted by the observer. Fig. (4.4-1) illustrates the output (y) 

I I I 

obtained from the gyroscope with sensing axis edong the e^ body axis. 

Figs. (4.4-2)-(4.4-5) give the angular velocity (Z^) and arngular 
displacement (Z 2 ) predicted by the fourth order observer model for the 
case when X^ = -0.5. The free response case is illustrated in Figs. (4.4-2) 
and (4.4-3), vhile Figs. (4.4-4) and (4.4-5) give the response vhen the 


wind acceleration is included. 



Figs. (4.4-6)-(4.4“13) illustrate the free response of the fifth 
order observer system with =-0.2 and - 0.5 respectively. Figs. (4.4-6) 
and (4.4-10) present the results for the angular velocity vtiile Figs. 
(4.4-7), (4.4-8), (4,4-11), and (4.4-12) present the results for the 
angular displacement. Figs. (4.4-9) and (4.4-13) contain the results 
for predicted bias (Z^) . 

Figs. (4. 4-14)- (4. 4-25) presait the fifth v rder observers response 
with wind acceleration for the case vhen X^ = -0.2, -0.5, and -0.7 
respectively. Figs. (4.4-14), (4.4-18), and (4.4-22) presoit the results 
for the angular velocity vhile Figs. (4.4-15) , (4.4-16) , (4.4-19) , 
(4.4-20), (4.4-23) and (4.4-24) present the angular displacement. The 
bias predicted from the fifth order daserver for the respective cases 
is presented in Figs. (4.4-17), (4.4-21), and (4.4-25). 
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Fig. (4.4-1) Actual Output For Gyroscope T'ounted Along X. Body Axis 
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Pig. (4.4-3^ Actual Results For Angular Displacement 
(4th Order Observer, u=0, X=-0._5) 
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(4.^i~12) Actual Results For Angular Displacement 
(5th Order Observer, u=0, A=-0,5) 
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4. 5 Discussion of Palloon Observer Results 

Figs. (4,4-2) and (4,4-4) indicate that the angular velocity 
results predicted by the fourth order observer differ significantly 
from those obtained from the gyroscope. The actual angular velocity 
of the balloon platform is different than that obtained from either the 
observer or gyroscope. Deficiencies in the ability of the gyroscope 
to reproduce the actual angular velocity are caused by mechanical 
inrensitlvity to sudden changes in angular velocity and bias. Errors 
in the angular velocity values predicted by the fourth order observer 
result from the form of the elements of the F and G matrices (see 
Sec. 4.3) and the magnitude of bias present in the output of the gyroscope. 

Fig. (4.4-3) shows that the angular displacement predicted by the 
fourth order observer model without wind input differs significantly 
from that obtained by integrating the output from the gyroscope (yi). 

The general trend of the angular displacement predicted by this model 
is sLmular to that obtained by integrating the output of the gyroscopes. 
However, the output of this model contains errors; e.g., the maximum 
difference between the two curves is 0,172°, This is significant in 
view of the fact that the maximum displacement predicted by the observer 
is 0.058° while that obtained by integrating the output from the gyro- 
scope is 0.23°. These errors can be attributed to the fact that: 

1. The effect of the wind acceleration is neglected, 

2, Bias is present in the plant output. 

Fig, (4.4-5) shows that the angular displaconents obtained from 
the fourth order observer including wind acceleration compare less 
favorably (with the Integrated gyroscope output) than those obtained 
from the previous model. The magnitude of the maximum angular displacement 
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prfidicted by this model is approximately 0,5°i while the majcimum 
deviation is of the order 0.6°. 

Figs. (4.4-6), (4.4-10), (4.4-14), (4,4-1^), and (4.4-22) show 
that the angular velocity results obtained from the fifth order observer 
model differ significajitly from the results (biased.) obtained from the 
gyroscope. This is true regardless of whether the effect of wind 
acceleration is included. However, as the magnitude of the repeated 
eigenvalue ()^j^) increases, the differences between the two curves decrease. 

Figs. (4.4-7)* (4,4-8), (4,4-11), and (4.4-12) show that the singular 
displacements predicted by the fifth order oteerver model which was 
developed byexcluding the affect of wind accleration, are in good 
agreement with those obtained by integrating the modified output from 
the gyroscope (yim) . The latter was obtained by subtracting the bias 
predicted by this model from the actual gyroscope output. These 
figures indicate that the difference between the two curves Increases 
with increasing time. Moreover, this difference decreases with in- 
creasing magnitude of The maximum displacement wlthXj^“-0,2 and 
-0.5 is of the order 0.05°; while the maximum error is 0.08° and 
0.02° respectively. 

Figs. (4.4-15), (4.4-16), (4.4-19), (4.4-20), (4.4-23), and (4.4-24) 
indicate that the angular displacements predicted by the fifth order 
observer which Includes the affect of wind acceleration are in good 
agreement with the values obtained by Integrating the modified output 
from the gyroscope. The difference between these two curves r«nains 
constant over the entire observed time period. Moreover, this difference 
decreases with increasing magnitudes of X^. '’■’he magnitude of the 

maximum error and displacement for all three cases is approxirvately 



0 , 2 5 *^, and 0.5*^ rospectlvely. 

Figs. (4J+--), (^'-.4-15), (4.4-16), (4.4-19), (4.4-20), (4.4-23), 
and (4.4-24) show that the angular displacements predicted by the fourth 
and fifth order observer models compare favorably regardless of the 
^^lue ofXj^. This Indicates that the predicted angular displacement 
is unchanged regardless which observer model is used. Thus, the 
fourth order observer model can be used to predict the angular displace- 
ment of the balloon platfom. 

?lgs. (4.4-3), (4.4-7), (4.4-9), (4.4-U), and (4.4-12) show 
that the angular displacements predicted by the observer models with- 
out wind input differ significantly from those results predicted by 
the models which include wind input, -This indicates that knowledge 
of wind accleration is necessary in order to obtain accurate results 
for the attitude of the platform. 

4.6 Conclusions 

This study has shown that, for a completely observable balloon 
systan, observer models can be constructed to accurately determine the 
angular displacement of the observational platform. Any errors in 
the predicted platform state are due mainly to errors in the balloon 
flight data (l.e, .acceleration and angular velocity data) as opposed 
to deficiency in the observer model. Although the results for the 
angular velocity are in poor agre«nent with those obtained from the 
gyroscopes, these deviations can be decreased considerably by proper 
choice of the eigenvalues. 

This study has also shown that the angular displacements predicted 



"by the observer models do not vary significantly with either the order 
of the model or the magnitude of the repeated elgen^lue. However, 
the results do vary significantly depending on whether the affect 


of wind acceleration is Included. 
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APPE24DIX A 


Balloon Translatioral Aoceleration Oonponenta 

In the 02186 of the lACAIE experiment, the balloon's position was 


tracked by radar; the tanslational ocnponents were obtained with res^>ect 
to the earth fixed axis shown in Fig. (A~l) . Ihe corresponding body 
axis for the balloon platform system is shown in Fig. (A-2) . Hie angle 
a meeisured between these two coordinate systans (Fig. A-2) is given as 
follows; i.e. , 


where 


a = 


t 

W^dt - tOs(X) 
o 


^ 0 ' 


(A-1) 


s spin ocnponent of amgular velocity obtained from gyroscope, 

= magnitude of earth spin, . . 

(7.2722 X lO"^ (rad * s"^) , 

X « latitude angle (0.5724 (rad) ) , emd 

~ initial value of at as measured by magnetometer. 

The velocity ocnfxinents of the balloon were obtained by numerically 

differentiating the (radar tracked) translaticxial oonponents. Hie velocity 

ocnponents of the ballocm (V , V ) measured along the balloon's body axis 

X y 

are given as follows; i.e. , 

Vi « - V^S(a) + Vy C (at), V 2 - C (a) + S (a) , (A-2) 

v4iere 


» balloon velocity component 
• 1 1 

eilong the 62 bo(^ axis, 

I I I 

» balloon velocity oonponent along the e^ body axis, 

V » ballocn velocity conponent along the e. earth fixed axis, 
Vy = balloon velocity oonponait along the earth fixed axis. 
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and a is as defined in Eq. (A**l) . 

Ihe balloon's translationed acceleration coi(X3nents alcxig the boc^ 

axis can be obtained by differ^tiating Eq. (A-2) with respect to time. 

Ihe resxilting equations are given as following; i.e . , 

“ Vj^S(a)+VyC(a)4t(-V^C(a)-VyS(a)) , (A-3) 

a, « V C(a)4V S(a)4«(-V S(a)4V C(a)), 

A jt y X y 

where 

a^ > translational acxeleration oonp o nent along the axis, 

a^ * translational acceleration coiponent along the e^' ' ' bod^ axis, 

• 

■ translational acx^eraticxi ocnponait along the earth fixed axis, 

V » translational acceleration ccnponent along the e. earth fixed axis, 

y ^ 

a * - QS(x) , and 

V , V and a are as defined previously. It Aould be noted that, in the 

* y 

above developoent, the earth's rotational effects are neglected. 
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Fortran Godina 


1 . Body Axis Accelerations 


1.000 c 

PROGRAM TO OBTAIN BOOT ACCEL. PROM 

2.000 C 

SEMsrrrviTY analysis 

■5.000 

ODTPOT ' INPUT N’ 

4.000 

INPUT N 

5.000 

DO 1 1=1. N 

6.000 

READ(103.2) T.P 

7.000 2 

FORMAT (2G) 

8.000 

READ(109.2)T.PD 

9.000 

READ(104.3) TO.XX.VX.AX 

10.000 3 

FORMAT (4G) 

11.000 

READ(107.5) T1.YY.VY.AY 

12.000 5 

FORMAT (4G) 

13.000 

A1 =AX*COS( P) +AY*SIN( P) ■»-(*VX*SrN (P) +VY*C03 (P) ) *PD 

14.000 

A2=AX*(-SIN(P) )4AY*C0S(P)+(-VX*C0S(P)-VY»SDKP) )*PD 

15.000 

WRITE(106.4) T.A1.A2 

16.000 4 

P0RMAT(3B14.6) 

17.000 1 

CONTINUE 

18.000 

STOP 

19.000 

END 
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2, Balloon Fourth Order Observer 


if. ^ 

OF POO,-^ QUALITY 


1.000 C 

2.000 
3.000 

4.000 

5.000 9 

6.000 

7.000 

8.000 
9.000 

’^.000 
1 : -000 
12.000 
13.000 

14.000 

15.000 
16.000 

17.000 

18.000 
19.000 
20.000 
21.000 
22.000 

23.000 

24.000 

25.000 

26.000 1 

27.000 

28.000 
29.000 
30.000 

31.000 

32.000 

33.000 

34.000 

35.000 4 

36.000 5 

37.000 

38.000 
39.000 2 
40.000 

41.000 

42.000 
43.000 

44.000 

45.000 

46.000 

47.000 

48.000 

49.000 10 

50.000 
51.000 


BALLOCBf 4TH ORDBR OBSTOYHl SYSTEM 
COrfOf/PCTT/A31 . A32, A41 , A42. 01 , 02, 03. 04, B32 
EXTSWAL PCT.OOTP 

DIMHJSIOM Y(5).DBRY(5).PRNT{5),AUX(16,8) 
OUTPOT’INPOT Y10.Y20' 

DATA (Y(I), I=1.4)/0.,0..0..0./ 

OOTPOT ’INPUT PrafP(I).I=1,4’ 

INPUT. (PRMT(I).I=1,4) 

OUTPUT 'INPUT EIO’ 

INPUT EIO 

A31 =-1.622 

A 32=1. 1926 

A41=6.1096 

A42=-8.1096 

B32=.0437 

04=-4*BI0 

03= (04*A51 -4*EIG**3 ) /A41 
G2=(-BI0»*4-A32*A41 +A42*A31 )/(-A32*A41 -»A42*A31 ) 
01 = (-6*EI0**2-A31 -A42+A42^»G2) / (-A41 ) 

OUTPUT 01.02.03.04 
NDIM=4 

OUTPUT'INPUT 1 TO RON' 

INPUT J 

IP (J .NE. 1) GO TO 9 
DO 1 1=1,4 
DERY(I)=0.25 

CALL HPCG(PRMT.Y,DERY,NDIM,IHLP,PCT.OUTP.AUX) 

STOP 

END 

SUBROUTINE PCT{X, Y.DERY, INO) 

DDQNSION Y{1).DERY(1) 

C0f«0N/PCTT/A31 , A32, A41 , A42. 01 , 02. 03. 04. B32 

C0MM0N/DAB/T1.THTD 

I?(INO .BQ. 0) GO TO 2 

FORMAT (20) 

FORMAT (30) 

READ(104,4)T1,THTD 

READ(103.5)T2,A1,A2 

DBRY(1 )=01*(-Y(4)+THTD)+Y(3) 

DERY ( 2 ) =02* ( -Y (4 ) +THTD) +Y (4 ) 

DERY 3)=G3*(-Y 4)+THTD)4A31*Y(1)+A32*Y(2)-B32*A2 
DBRY(4) =04* v-Y(4) +THTD) +A41 *Y ( 1 ) +A42*Y(2) 

RETURN 

EiHD 

SUBROUTINE OUTP(X.Y,DERY,IHLP.NDIM,PRMT) 

C0W0N/DAB/T1.THTD 

DMaiSION Y(1),DBRY(1),PHMT(1) 

WRITE( 1 06, 1 0)X. Y(2) , Y(4) . THTD 
P0RMAT(4E13.5) 

RETURN 

QID 
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PASS !S 

OF POOR QUALITY 

3. Balloon Fifth Order Observer 


1.000 C 

2.000 
3.000 
4.000 
5.000 9 
6.000 

7.000 

8.000 
9.000 

10.000 

11.000 

12.000 

13.000 

14.000 

15.000 
16.000 

17.000 

18.000 

19.000 

20.000 
21.000 
22.000 
23.000 

24.000 

25.000 
26.000 
27.000 1 
28.000 
29.000 
30.000 
31.000 

32.000 

33.000 

34.000 

35.000 

36.000 4 

37.000 5 

38.000 

39.000 

40.000 2 

41.000 
42.000 

43.000 

44.000 

45.000 
46.000 

47.000 

48.000 

49.000 

50.000 
51.000 10 
52.000 
53.000 


PR0®AM TO TEST 5TH ORDER OBSHIVSI SYST01 WITH INPUT BIAS 
C0W0N/PCTT/A31 , A32. A41 , A42, G1 . G2, G3, G4, G5, B32 
EXTERNAL PCT,OUTP 

DIMiNSION Y(5),DERY(5),PRMT(5),AUX(16.8) 

OUTPUT 'INPUT Y10,Y20,Y30' 

DATA (Y(I), I=1,5)/0.,0.,0..0.,0./ 

OUTPUT 'INPUT PRMT(I).I=1,4’ 

INPOT, (PHMT(I),I=1,4) 

OUTPOT 'INPUT EIG' 

INPUT EIG 

A31 ^1.622 

A32=1.1926 

A41=8.1096 

A42=-8.1096 

B32=.0437 

G2=(5*EIG**4-A42*A31 4A41 *A32) /(-A42*A31 +A41 »A32) 

G1 = { 1 0*EIG**2-A42*G24A31 4A42 ) /A41 
G5^BIG**5/ ( A42*A31 -A41 *A32 ) 

G4=-5*BIG-G5 

G3=(-1 0*EIG**3+A31 ^KJ4-^(A31 ->A42) ^5) /A41 

OUTPUT G1,G2,G3.G4,G5 

NDIM=5 

OUTPUT 'INPUT 1 TO RUN' 

INPOT J 

IP (J .NE. 1) GO to 9 
DO 1 1=1,5 
DERY(I)=0.2 

CALL HPCG(PRMT,Y,DBRY,NDIM,IHLP,PCT,OUTP,AUX) 

STOP 

END 

SUBROITPINE PCT (X. Y, DERY, INO) 

DIMMSION Y(1),DERY(1) 

C0fW0N/PCTT/A31 , A32, A41 , A42, G1 , G2, G3, G4, G5, B32 

C0TW0N/DAB/T1.THTD 

IP(INO .BQ. 0) GO TO 2 

FORMAT (2G) 

FORMAT (3G) 

READ(104,4)T1,THTD 

READ(103,5)T2,A1,A2 

DERYd )=Y(3)-tG1*(-Y(4)-Y(5)+THTD) 

DERY ( 2 ) =Y ( 4 ) -K52* ( - Y ( 4 ) - Y ( 5 ) +THTD ) 

DERY(3)=A31 *Y(1 )+A32*Y(2)4G3*(-Y(4)-Y(5)+THTD)-B32*A2 
DERY ( 4 ) =A4 1 * Y ( 1) 4A 42* Y ( 2 ) 404* ( - Y ( 4 ) -Y ( 5 ) +THTD ) 
DERY(5)=G5*(-Y(4)-Y(5)4-THTD) 

RETURN 

SUBROUTINE OUTP(X,Y,DERY, IHLP.NDIM.PRMT) 

C0MM0N/DAB/T1,THTD 

DIMENSION Y(1), DERYd ),PRMT(1) 

WRITBd 06, 1 0)X, Y(2) , Y(4) , Y(5) . THTD 
PORMAT(5E13.5) 

RETURN 

END 
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4. Hamming* Predictor Corrector 


1.000 c 

?.ccc c 

J.CCC PrPFOUTINE HPCG(PPPT,y.rFFY,NriM.IHlF.FCT.OUTP,AUX) 

4.0C0 C 
F.OCO C 

6.000 PJMFNFIOK PFPT( 1 ) . Y( 1 ) . DFRY ( 1 ) . AUX (1 6, 1 ) 

7.000 N=1 

e.OOO IFIF^O 

9.000 XrFPKTfl) 

10.000 H=PPPT(Y) 

11.000 FFPT(5)*0 

12.000 rO 1 1=1 NDIM 

13.000 Auxne. i)=o. 

14.000 Anc(15. 1)*DFPY(I) 

15.000 1 AUXn.I)=Y(I) 

16.000 IF(H*(PPFT(2)-X))3,2,4 

17.000 C 

18.000 C FFFCR RETURN? 

19.000 2 IH1F=12 

20.000 GOTO 4 

21.000 3 IH1F=13 

22.000 C 

23.000 C COMPUTATION OF TERY FOR STARTING VALUES 

24.000 4 IK0=1 

25.000- CALL FCT(X,Y,rEFY,INO) 

26.000 C 

27.000 C PECOPDING OF STARTING VALUES 

28.000 CALL OtTP(X, Y.TERY, IHLF,NDIM, PRMT) 

29.000 IF(PRMT(5))6,‘^,e 

’0.000 ^ IF(IHIF)7.7,6 

’1.000 6 FFTITN 

’2.000 7 ro 8 I=1,NriM 

’3.000 8 Ara(8.!)=rFPY(I) 

34.000 C 

’5.000 C COMPUTATION OF AUXfP.I) 

36.C0C IFV=1 

37.000 GOTO 100 

38.000 C 

39.000 9 X=X4fi 

40.000 ro 10 i=i.NriM 

41 .000 10 AUX(2 I)*Y(I) 

42.000 C 

43.000 C INCPFMFKT H IS TESTEP PY MEANS OF BISECTION 

44.000 11 IHLF=IHLF4l 

45.000 X*X-H 

46.000 PO 12 I=1,NPIM 

47.000 12 AUX(4,])=AUX(2,I) 

48.000 H=H 

49.GOC N=1 

50.000 ISV*2 

51.000 GOTO 100 

52.000 C 



ORIGirj.U PAGE fS 
OF POOR QUALJr/ 


108 


5?.OCO 
8/3 CCC 
f?.CCC 

^e.occ 

^7. CCC 

coc 

5C.CCC 

60 . CCC 

61 . CCC 

62.000 C 
6J.0CC C 
64. COO 
65 . OOC 

66.000 

67.000 

68.000 
69.000 
70. COC C 

71 . 000 C 
72.0C0 
7‘* . OCO 
74 . OOC 

75. CC0 C 

76. CCC C 
77 . OOC 

76.000 
"^o.COC 
80. COO 

81 .000 
82 . OCC 
8; . OCC 
84 . OOC 
85-000 
86. CCC C 
87 . OCC 
88.000 
89.000 

90 . CCC 

91 . COC 
90.000 
91 . 000 
94 . COO 
95. OOC 
96 . CCC 
97.000 
98. OCC 
QO. 000 

1 00 . OCC 
101 .000 
102.000 
105.00c 

1C4 .OCC 
1C5-CCC 
1C6.0CC 


1 5 


1 4 


X=X4H 


TFC^C 

CALI FCT(x.y.rFFy ikc) 
K^2 

rc 14 1^1 KriF 
Arx(2. I )=Y(I ) 

AUX(9 I)=rFFY(I) 

GCTC ICO 


CCFFUTATICK CF TEST VAICE FELT 
15 FELT^O. 

FC 16 I^I.NFIM 

1 6 FF1T=FE1T4AEX ( 15 , I )^AE£( Y(I )-AUX (4 , I ) ) 

FFLT=.C6666667*FEIT 
GO TC 19 

17 TF(IHLF-10)1 1 . 18, 18 

NO 8ATT8FACTCFY ACCURACY AFTER 10 PI8FCTJCNS. ERROR MESSAGE. 

18 THIF=11 
X^X4H 
GOTO 4 


'"FFFF IS SATISFACTORY ACCURACY AFTFF lESS THAN 11 PISECTICNS 
10 X=X 4 H 
INC=0 

CALL FCT(X, Y, FFFY. INC) 

FC 20 1=1 . NFIM 
AUX(5. I)=Yn) 

20 A.n((1C. I)=FFFY(I) 

N=5 

ISV=4 
GCTC 100 

21 N = 1 
X=X4H 

INC=0 

CAll FCT(X. Y.FEFY. INC) 

X=PFMT(1 ) 

FO 22 1=1 ,KFIM 
A UX (11 I ) =FER Y ( I ) 

220Y(I)=AUX(1 , I)4H*( .375*AUX(8, I)4.7916667*AUX(9, I) 
1-.2083553*AUX(10. I )4.041 66667^^FERY(I ) ) 

23 X=X4H 
N=N4l 

INC = 1 

CALI FCT(X, Y,FEFY, INC) 

CAII OUTF(X, Y. FFFY, IHIF,NFIM, FFMT) 

IF(FFMT(5) )6, 24, 6 

24 TF(N-4)25, 200, 2CC 

25 FC 26 1=1 . NFIM 
AUXiN, I)=Y(I) 

26 AUX(N47. I)=FFFY(I) 

IF(N-3)27. 29 , 2CC 


1C9 


1C7.CCC 

Hb.CCC 

1CO.CCC 

1 10. CCC 

111 . ccc 
1 ir.CCC 
1 1 J.OOC 

114. ccc 
1 1 *^.ccc 

1 ’ 6.CCC 

1 1 7 . ccc 
1 1 e.ccc 
1 1 ccc 
1 2 C . ccc 
1 ri .ccc 
1 2?. ccc 
12 ?. ccc 
1 24 .ccc 
1 2 *^. ccc 
i 2 e.ccc 
127. ccc 
12b.00C 
129.0CC 
1?C.CCC 

1 ■* 1 . ccc 
1 : 2 . ccc 
1 ?: .ooc 

1 ;4.CCC 

1 ??.ccc 
1 ’t. coc 
1:7. ccc 
1 rc.ccc 

1 39 . ccc 

1 4C.CCC 

141. ccc 

142. ccc 
1 4 ?. ccc 
1 44 . ccc 
1 4 . ccc 
1 4 e . ccc 
1 4"-. ccc 
148. ccc 

1 40. ccc 

1 9C.CCC 

1 91 .ccc 
1 ^ 2 . ccc 
1 9 ? . ccc 
1 94 . ccc 
1 '^ 9 . ccc 

196 -CCC 
1 97. ccc 
1 96.CCO 
1 99. ccc 
1 tc.ccc 


ORIGIN.^L TAG.* JS 
Q OF POOR QUAL.'Vr 

27 :>C 28 Ir1.K^]^ 
rFLTrAl’X(o. I)-tAl'X(9. ] ) 
rFlTrCFlT-fTFIT 

28 y {] )^Ai'X ( 1 . ] )-f.????:?:*p*(Ai’> ( 8 . ! ) 4 riLT-fAa(ic, ] ) ) 

GCTC 2? 

c 

20 PC ?C 1^1 . KHP 

rFiTrAuy(9. i)4Arx(ic, i) 

^FITr^FIT-^^FIT^^FIT 

?c y ( I ):rAry ( 1 . i)4.?79*hmai’X( 8. ] )4rFiT4Ary ( 1 1 . i ) ) 

GCTC 27 
C 

C THF FCIICWING FAFT CF Pl’FFCU'I'JKF FFCG CCNFUTF8 Py CF 

C FITGF-KTTTA FFTPCF PTAFTING VAIPFP FCF TPF POT PFIF- PTAFT] KG 
C FFFr’CTCF CCFFFCTCF PFTPCr. 

ICC rc 1 C 1 1 ^ 1 ,NriF 

z^p^Apy (N 4 ^. ] ) 

Apy('^. ])=z 

1C1 y(i)=Apy(N i)-t.4*L 
C Z IP AN AUXIllAFy PTCEAGF ICCATICK 

c 

ZrX4 4*P 
INC^C 

CAII FCIlZ. y, TEF^, INC) 

rc 1 C 2 iri.KriF 

Z=H*rFEMI^ 

Al’X lb, I )=Z 

1C2 yiI)=AUX(N. 1)4.2969776*AUX(9. I ) -f . 1 967596*Z 
C 

Z^y-f. 4997772 *P 

INC=C 

CAiL FCT(z. y. rFpy, ikc) 

rc 1C7 1 ^ 1 , Nr IK 

ZrP*rFFy(i'» 

AlT(( 7 .I)^Z 

IC"* Yn)rAPy (N. I)4.2181CC4*A17(9. I)-'^.C9C969*APy (6. 1)47.872869*2 
C 

Z^X4P 

IKC^C 

CAII FCTfz. y. rFpy, INC) 

rc 1C4 iri.NriK 

iC4cy(i)^Arx(N. 1 ) 4 . i74':6C7*A.ry (9. i)-.57i^ec7*Ary(6 i) 

Ml . 2 C 9 ‘=‘* 6 *APy ( 7 . iw. 1 71 1 e 48 *H*rEFy( I ) 

GCTC (0.17 19 .ri).IFW 

c 

c FCPPIFIF FFFAK-FCINT FCF LINKAGE 

C 

C PTAFTING VALPEP AFF CCKIPTEF. 

C NCV PTAFI HAKKINGP K.CriFIEr FREri CTCE- CCFFECTCF KFTKCr. 

2CC IPTFFr* 

2C1 IF(N-8)2C4,2C2, 2C4 
C 

C N^8 CAPPEP TPF FCWP CF APX TC CHANGE THFIF PTCFAGE LCCATICN^ 










IIU 


f 


161 .cor 

16?.CCC 
16?.CCC 
1 64 .OCC 
1 6?.CCC 
1 66.COC 
167.000 
160.000 
1 69.000 

170.000 

171 .000 
172.000 

173.000 

174.000 

175.000 

176.000 

177.000 

176.000 

179.000 

160.000 
161.000 
162.000 
163.000 

1 64 . 000 

165.000 
166.000 

167.000 
1 66 . 000 
16<^.C00 
1 90 . 000 

1 01 . 000 
lo^-.OOC 
1 93 . 000 
1 Q4 . OCC 

1 OP . 000 
1 06 . 000 
1 97 . 000 
1 96 . 000 

1 99.000 

200 . 000 
201 .000 
202.000 

203 . 000 
204 . COO 

205.000 

206.000 

207.000 
2CS.CC0 
209. OCC 
210. OCC 

21 1 .000 
21 2. OCC 
21 3. OCC 
214.000 



21 

0 

0 

C 

0 


ORlG^'i'V- PAGH iS » 

OF POOR QUALITY 

202 ro 203 N^2.7 
ro 203 i^i.Nriy 

Al’>’ (N-1 .])=AVy(U, I) 

203 Ary(N46, I)=ArX(N.f7. I) 

Nr7 

K ir?P THAK 6 CAl'FF? FN 1 TO GIT H 

204 N^N4l 

CCFFUTATICN OF KFXT VFCTCF Y 
DC 205 I=1.fJDIN, 

ArX(N-1 . ])=Y(I ) 

205 ADX(K46, I )=rFFY(I ) 

X=XhH 

206 ISTEFrlSTFF^I 
DC 207 I = 1,NDIFi 

CDIIT=AUX(N-4, 1 )4l . 333333*H*(AUX(N46, I )4AI’X (N-ft , I )-AUX (F45, 1 ) 
1AUX(N44, l)4ArX(N44, I)) 
y (I )=DFIT-. 92561 96*AUX (1 6, I ) 

207 A rx ( 1 6 , J ) ^DFIT 

FFFDICTCP I? NOW GFKFFATFD IN FCW 16 OF AUX, l^CDIFIED FFFDICj 
1? GFNFFATFD IN Y. DELT NFAN6 AN Al’XIIIAFY FTCFAGF. 

INC^O 

CAIl FCT(X, Y, DEFY, INC) 

DFPTVATIVF OF NCDTFIFD FFEDTOTOF IS GFNFFATFD IN DEFY 
DC 206 I^I.NDIK 

ODFIT^. 125*(Q.*AUX(N-1 , I)-AUX(N-3. I ) 43 . (DEFY f I )4A I’X (N4 6 . I) 
1Al7(N.t6, I)-ArX(N49. I)) ) 

AUX(1 6, D^Al’Xd 6. I)-DFIT 
206 y (I )rDFIT4. 074 3601 7'*ArX( 1 6. I ) 

TEST WHFTHFF H NTST FF KAIVFD OF DCUFIFD 
DF1T=C. 

DC 209* 1^1 .NDIW 

209 DF1T = DFLT4AIX(1 5. I ) ♦AES ( A FX ( 1 6 . I ) ) 

GO TO 210 

H NUST NOT FF HAIVED. THAT WFANS Y(I) AFF GCCD. 

210 INC^I 

CALL FCT(X, Y.DEFY, INC) 

CALL 0UTI{X,y,DEFY, IHLF, NDIM , FF^T ) 

IF(FFNT(5))212,21 1 , 21 2 

211 IF(IHLF-11 )213,212,212 

212 FETUPN 

213 IF(E*(X-PFMT(2)))214.212.212 

214 IF(AFS(X-FPNT(2) )-. 1*AFS(H) )21 2, 21 5. 21 5 
5 GO TO 201 


H COULD FF DOUFLED IF ALL NFCESSAFY FFFCFFDING VALUFS AFF 
A VA ILAFTF 

216 IFfIFIF)20i .201 ,21'7 


?1?.000 

?16.CCC 

?17.C0C 

pie.ooc 

?19-OOC 
P20.0CC 
. OCO 

???.ooo 

22? . COC 
224. OCC 

225.000 

226.000 
227 • 000 
228.000 

229.000 

230.000 

231 .000 

232.000 

233.000 

234 . 000 

235.000 
23f .000 

237.000 

238.000 

230.000 

240.000 

241.000 
24?. OCO 

243.000 

244.000 

245.000 

246.000 
247 -000 

248.000 

249.000 

250. 000 

251 .000 
2 ^ 2.000 
253 000 

254.000 
255-000 

256.000 

257.000 

258. 000 

259.000 

260.000 
261 .000 
262,000 
263.000 

264.000 

265.000 
266.000 
267.000 


original Pag?" ■ 

217 IF(F-7)201 . 218. 218 QUAln Y 

218 TFaPTFP-4)20^219.219 
?lo TFOrsTF'^FP/? 

TFfTFTFF- TFOO IF0r)?0i , 220. 201 
220 H=-F 

IH1F=-IHLF-1 

I?TFF=0 

ro 221 iTi.FnM 

AUF(N-1 . 1)rAra(K 2,1) 

AUX(N-2.I)*^ArF(N 4,1) 

AUXfN-3-I)=AUX(N-6,I) 

ArafK46. I)sAI’X(N45. I) 

AUX(N45 - I)=AU2(»43. I) 

AUX(N44.I)*AUX(N4l,l) 

DELTsArX(N46.I)4AUX(N45, I) 

DEIT s-DEIT 4 DEIT 4PE1T 

2210 AUX( 16 . 1)*8.962963*(X(I)-AUX(N-3,I))-5.561111*H*(r?ERX(I)4l'ElT 
14AUX(H44. I)) 

GOTO 201 


F MOST PF HAIVEr 

222 TFIF=THIF4l 
IF(TF1F-10)2??,223,210 

223 P^F 
TFTFP^O 

ro 224 I=>1.NriF 

0y(T)=.003?06?5*{80.*AUX(N-1,I)4l35.*AUX(N-2,I)440.*AUX(N-3, 1)4 
1AUX(N-4.T))-.117l875*fAUX(NH6.T).-6.*ArX(N45. 1)-AUX(N44. 1 ) )*H 
OA FXCF-4, 1 )*. 00300625* (1 2. *A UX(N- 1. 1)4l 35- *Ara (P-2, 1)4 
1108*AUX(P- 3, l)->Ary(P-4, 1)) . 0234375* (AUX(P46. I )4l 8- *AUX(N45. I ) 
29. *AUX(N44. I))«F 
ArX(N-3- I)=AlT7(P-2. 1) 

224 AUX(N44. T)rArafN45.I) 

x=-y-H 

PFITrX- (E4F) 

INOs^O 

CALI FCT (PELT. Y, DEFY INO) 

DO 225 I=1.NriF 

AUX(N-2. T)=y(I ) 

AlTX(N45.I)=rEPY(]) 

225 Y(I)=AUX(K-4,I) 

PElTtrLFlT- (H4H) 

INOeO 

CALL FCT (PELT, y, PER Y, INO) 

PC 226 I*1,NPIM 
PELT=A0X{N45,I)4AFX(N44,I) 

PELT^PELT4PELT4PELT 

0ArX(16,I)s8.96?S63*(Anr(N-1,l)-y{l))-3.361111*H*(AUX(K46.I)4PEL 
UPFRYd)) 

226 AUX(N43, I)=PFFY(I) 

GOTO 206 
FPr 
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